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It is shown that anisotropic spin chains with gapped bulk excitations and magnetically ordered 
ground states offer a promising platform for quantum computation, which bridges the conventional 
single-spin-based qubit concept with recently developed topological Majorana-based proposals. We 
show how to realize the single-qubit Hadamard, phase, and ty/8 gates as well as the two-qubit 
CNOT gate, which together form a fault-tolerant universal set of quantum gates. The gates are 
implemented by judiciously controlling Ising exchange and magnetic fields along a network of spin 
chains, with each individual qubit furnished by a spin-chain segment. A subset of single-qubit 
operations is geometric in nature, relying on control of anisotropy of spin interactions rather than 
their strength. We contrast topological aspects of the anisotropic spin-chain networks to those of 
p-wave superconducting wires discussed in the literature. 
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I. INTRODUCTION 

Quantum control of electronic systems is usually re- 
stricted to the realm of very small structures, wherein 
qubits are engendered by elementary quantum degrees of 
freedom like electron spin (or, more generally, Kramers 
doublet). Notable exceptions are provided by macro- 
scopically coherent superconducting qubitJ^ and anyons 
with non-Abelian braiding statistics!^ Such "nonabe- 
lions" may arise in certain correlated topological states 
of matter 3 (most prominently, the v = 5/2 fractional 
quantum Hall effect) or be fabricated in topological- 
insulator based systems 4 and even rather conventional 
semiconductor-based heterostructuresP The latter re- 
quire tailoring an appropriate combination of spin-orbit 
interactions, magnetism (or magnetic field), and su- 
perconductivity in order to furnish a topological state 
with non-Abelian Majorana fermions. The localized 
and spatially separated Majorana fermions can endow 
the ground state with necessary degeneracy 6 and robust 
quantum coherence^ which in turn provide a basis for 
topologically-protected quantum computationP 

Despite a promise for solid-state quantum computa- 
tion, both the superconductor-based and topological (al- 
most always Majorana-based) platforms are not without 
serious challenges, however. The former suffers from rela- 
tively fast (while poorly understood) decoherence, which 
is perhaps related to its macroscopic nature and non- 
superfluid parasitic low-lying degrees of freedom, while 
the latter has in general only a subset of quantum gates 
that is topologically protected, which needs to be sup- 
plemented by more traditional schemes (for example, 
based on interaction between Majoranasp' coupling to 
flux qubits, 10 or an rf measurementPH) , in order to real- 
ize a universal set of gatesP^ Overcoming these obstacles 
poses tremendous technological challenges. 

At the same time, steady progress is being made in 
harnessing single-electron spins in semiconductor quan- 
tum dots as natural building blocks for universal quan- 
tum computation.^ While a serious problem is posed by 



nuclear-spin induced qubit decoherence, record- long co- 
herence times have been achieved, 14 which are four orders 
of magnitude larger than those in early experiments.^ 

In this paper, we are proposing a hybrid platform 
for universal quantum computation, which is based on 
anisotropic spin chains. A chain can consist of a finite 
sequence of several appropriately engineered and tun- 
able quantum dots or a section of magnetically ordered 
one-dimensional wire. In the absence of magnetic fields, 
the degeneracy of such a spin chain is guaranteed by 
time-reversal symmetry. In the presence of weak mag- 
netic fields, the degeneracy can survive up to a critical 
field, if the spin chain is sufficiently long on the scale 
of an appropriate coherence length. A nonlocal topo- 
logical character of spin-chain qubits may be inferred 
from the sensitivity of the quantum coherence within 
the ground-state subspace to the boundary conditions. 
Furthermore, a formal similarity of the nonlocal spin 
qubits in open chains with Majorana end states in topo- 
logical superconducting wires is revealed by a Jordan- 
Wigner transformation,^ allowing one to tap into the 
existing pool of ideas on Majorana-based quantum com- 
putation (while appreciating fundamental differences be- 
tween spin- and superconductor-based networks of such 
wires, as discussed in the Appendix). 

We will construct geometric gates on our spin systems 
that mimic braiding of nonabelions, while the readout 
can be performed via standard spin-to-charge conversion 
schemes^ on individual sites or dots. The readout can 
also be accomplished by fusing the spin-chain qubits and 
measuring spin torque, in analogy to the Josephson cur- 
rent readout of Ref. El While drawing inspiration from 
both worlds of (spin-based and topological) quantum 
computation, we are motivated by a versatility for practi- 
cal realizations of such hybrid systems as much as bridg- 
ing and synthesizing diverse theoretical ideas that have 
been recently developed in different subareas of quantum 
computation. 



II. SINGLE-QUBIT GATES 

The nonlocal spin qubits arising in anisotropic spin 
chains are defined (along with a review of necessary in- 
troductory material) in the Appendix. As a key result 
of this section, we will construct a phase gate on such 
qubits, which is analogous to braiding of two vortices in 
chiral p-wave superconductors^ or two Majoranas in su- 
perconducting nanowire networks P While our phase gate 
will be geometric in spirit, no braiding is necessary, how- 
ever. In fact, braiding of our Majoranas as in Ref. |8] 
would not generate a nontrivial gateP^ (See also Ref. [18] 
for braiding- less gates on Majoranas in superconducting 
wires coupled to single-electron quantum dots.) 

Before proceeding, let us specify our single-spin basis: 
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where (0,0) are polar angles (over the Bloch sphere), 
which define spin direction of a coherent-spin state, and 
It (I)) = <W(I) are t ne usual spin- 1/2 up or down eigen- 
state along the z axis. The singularity (Dirac string) of 
our basis is along the south pole on the Bloch sphere. 
The corresponding Berry phase^ is given by 
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In this section, we will construct single-qubit gates for 
an Ising chain. To provide an explicit demonstration, we 
start with a three-site segment [which is readily gener- 
alizable to longer chains, as sketched in Fig. [2^b) below; 
in practice, however, we would be interested in a meso- 
scopic situation, such that the chains are long enough on 
the scale of the appropriate correlation length £, to have 
the degenerate ground state, but not too long to avoid 
inevitable decoherence effects and errors associated with 
the control of magnetic anisotropics] described by the 
following time-dependent Hamiltonian: 
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When m = x and hi = 0, the chain is in some superpo- 
sition of the ground states |=>) and |^=), Eq. ( |A8| ). 

We first discuss a straightforward way to initialize the 
system in a Greenberg er-Hor ne-Zeilinger (GHZ) state 
|+), where [see also Eq. (|A12|)] 
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are the logic states for three spins (and analogously for 
N spins) . Ramping the uniform field hi —> h in the chain 
with m = x up to a large positive value establishes the 
unique ground state |ff). Slowly then turning off h\(t) 
to zero leaves the first site in the \\) = + \<—))/y/2 
state, while the second and third sites remain pinned 
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FIG. 1. (color online). Initialization of an iV-site (N — 3 as 
drawn) chain in a Greenberger-Horne-Zeilinger state |+). See 
text for details. 



in the same up state by the magnetic field along the z 
axis. Assuming, that h is sufficiently large, we can ne- 
glect quantum fluctuations of the pinned spins, so that, 
since (af ) = 0, the first site is completely decoupled from 
the rest of the chain. Turning then h2(t) slowly off, the 
second site picks coherently the orientation of the first 
site along the x axis, while the last site still remains de- 
coupled by the pinning field. This gives the following 
state for the first two spins: (!—>•—>•) + \<—<—))/y/2. In 
order to see this explicitly, we write the initial state for 
the first two spins as 



(5) 



where \s) describes the state of the second spin (effec- 
tively decoupled from the third spin). Since the interac- 
tion between the two spins is of the x-Ising form, which 
cannot flip the first spin along the x axis, their time- 
dependent state will have the form 
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with the initial condition I s+(t)) = |s_ (t)) = \s). Di- 
rect substitution of Eq. (|6| into the Schrodinger equa- 
tion then immediately shows that the second spin \s±(t)) 
evolves adiabatically according to the effective magnetic 
field h e ff = ±Jx + ^2 (^) z ? thus ending up in the entan- 
gled V% state. Note that in this process, 
Berry phases Q vanish in our basis, since 5(f) = 0. 

Finally, turning off hs(t)^ we establish the GHZ ground 
state |+), Eq. Q. The schematic of this procedure, 
which "zips" the state into |+) by gradually turning 
off magnetic field on sites 1, 2, etc., is shown in Fig. [l] 
Any initial quantum coherence in the logic space would 
be wiped out by this initialization. It is, furthermore, 
clear that this initialization works for an arbitrary num- 
ber of spins in the chain J 2 ^ 

We now implement the phase gate 
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on this three-site system as follows [see Fig. [2|a) for a 
schematic!. Let us start in the |=>) state with m = x and 
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hi = in Eq. We first ramp field hi up to a large 
positive value, which tilts the first spin in the xz plane 
from the x toward the z direction, contributing a vanish- 
ing Berry phase Q as 8(j) — 0. Then m is adiabatically 
rotated in the xy plane from x to y. This step rotates the 
remaining two spins in the xy plane from the x direction 
to the y direction, which contributes the Berry phase 7r/4 
per spin, since cos# = and S(j) = tt/2. The third spin is 
then biased by a large positive field h% (no Berry phase). 
The second spin is now effectively decoupled and in the 
state 
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When the field hi is slowly turned off at the next step, 
we thus arrive at the state 
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for the first two spins. Rotating m back to x and turning 
off hs returns us to the original Hamiltonian, while the 
wave function transforms to the state 
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neglecting the aforementioned Berry phase, as well as the 
dynamic phase, which are the same for the other initial 
state |^=). Repeating these steps for the |^=) initial state, 
we obtain the transformation: 
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Putting Eqs. (101) and (11) together, we get precisely the 
phase gate (|7) in the logic-space basis Q. The exten- 
sion of this scheme to longer chains is straightforward, 
which is crudely sketched in Fig.[2|b). We note in passing 
that this geometric procedure is analogous to Majorana 
braiding in superconducting wire T junctions,^ which ef- 
fectively flips the superconducting phase by tt. In our 
spin chain, rotating Ising spins by tt/2 in the xy plane is 
in fact fully analogous to changing the superconducting 
phase by 7r for the Jordan- Wigner fermions (see Sec. IV). 

Partial rotations of the Ising interaction axis from x 
to some direction m in the xy plane, at the intermediate 
step of the above procedure, could produce an arbitrary 
desired rotation around the z axis on the Bloch sphere of 
the logic space. Namely, for m = (cos0, shi(/>, 0), we find 
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The case <j) = 7r/4, in particular, realizes the ir/8 gate, 
which in the proposal of Ref. [TT] would require bring- 
ing Major anas together and lifting the degeneracy in the 
logic-space basis. 

Finally, a Hadamard-like gate (or likewise an arbitrary 
rotation around the x axis on the Bloch sphere of the 
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FIG. 2. (a) Application of the phase gate in the |=b) basis 
to a three-site spin chain \=>), as explained in the text, (b) 
Crude schematics for a longer chain: The blue regions show 
the "topological" Ising sections in the absence of magnetic 
field, with degenerate ground states. The red regions are or- 
dinary nondegenerate spin-polarized sections subjected to a 
large magnetic field in the z direction. The white arrows show 
the direction of the Ising interaction axis at each stage. The 
topological region encoding a qubit is first pushed to the right 
by raising the magnetic field to its left and lowering it to its 
right. The key operation is rotating the Ising axis from the 
x direction to y in the second step. The Hamiltonian returns 
to its original form in the last step, while the qubit undergoes 
the phase-gate transformation Note that there need not 
be a full control of individual sites, only anisotropics over the 
relevant segments of the entire chain. 



logic space) is easily realized by applying a magnetic field 
pulse (p = JhAt along the x axis to, say, the first spin of 
the chain, so that 
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For if = 7r/4, this gives 
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The Hadamard gate can be then obtained from this with 
the help of two phase gates: 
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Notice a certain symmetry between the bases (i) 
an d (h) {|=>), In the basis (i), we 

can implement an arbitrary phase gate (12) by rotat- 



ing the Ising interaction axis m, as described above, and 
a Hadamar-like gate (14) by the magnetic field pulse 



ip = 7r/4. In basis (ii), an arbitrary phase gate (13) is 
realized by a magnetic-field pulse, while the Hadamar- 
like gate, Eqs. (10) and ( pT| , by the tt/2 rotation of the 
Ising interaction axis in the xy plane. According to the 
preceding discussion, we can thus realize the desired set 
of the pure Hadamar, phase, and tt/S gates in either of 
the two bases. 
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III. TWO-QUBIT GATES 



For the two-qubit gates, consider two Ising spin c hains 
put in parallel, each described by Hamiltonian (A14) 
with h = 0. These qubits can be implemented as mo- 
bile topological sections embedded in longer spin chains 
and separated by nontopological regions as sketched in 
the example of Fig.[2|b). We can implement a two-qubit 
operation by pulsing an x-Ising exchange interaction J 
between, say, the first sites in each qubit. In the basis 
{|=>*), |^=)} for each qubit, the transformation amounts 
to 
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where 



JAt. For the special case of 
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Eqs. (16) result in the following gate in the logic-space 
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At this point, it is instructive to make a brief digression 
and recall the braid-group representation of four vortices 
in chiral p-wave superconductors P The vortices define 
four Majoranas, {71,72,73,74}, which can be combined 
into two complex fermions c\ = (71 + 272)/ '2 and C2 = 
(73 + 274 )/2. Braiding on a plane is accomplished by 
three generators, {Ti 2 , T23, T 34 }, respectively exchanging 
vortices 1 and 2, 2 and 3, and 3 and 4. These generators 
are represented by 



Tij = t(T 13 ) = e^/ 4 = 
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which introduces a relative sign between two Majoranas 
upon their exchange: 



li 
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such that ji — ^ 7j and jj — 7^ upon braiding vortices 
i 7^ j. 21 Written in the basis {1, cj, cj, cjc^HO), the braid 
group representation (18) is given explicitly by 6 
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States {l ? c{c2}|0) and {cj,c4}|0) have different fermion 
parities and are thus not mixed by braiding, which con- 
serves the parity. 6 We notice that T23 coincides with 
transformation (17) in our two-qubit logic space, while 
T12 and T34 simply reproduce phase gate ^ performed 
individually on one of the two spin chains. 



Since such standard braiding operations (20) are not 
sufficient for a universal quantum computation, we 
switch the logic-space basis from {|+), |— )} to {|=>), \<=)} 
and supplement transformation ( [l6| ) (0 = — 7r/4) with 
single-qubit phase gate (13) (ip = 7r/4) for each qubit. 



This gives (up to an overall phase) 
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The CNOT gate together with the Hadamard, phase, 
and 7r/8 gates provide a fault-tolerant universal set for 



quantum computation. 
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IV. FUSION AND READ-OUT 
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which is the controlled-Z gate. Sandwiching it by the 
Hadamard gates applied to the first qubit realizes the 



Bringing two topological segments in contact along a 
spin chain leads to interaction between the end Majo- 
ranas. Let us model it by the following Hamiltonian: 
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Here, two anisotropic sections i = 1 . . . M and i = M + 
1 ... AT have Ising interactions along m = (cos 0, sin (/>, 0) 
and x, respectively (recall a™ ~ &i ' m ) •> interacting with 
an isotropic (real- valued) XX coupling t in the xy plane 
at the junction, i = {M. M + 1}. When t -)> 0, the two 
sections are disconnected, each producing two Majorana 
sites at the ends, as explained in the Appendix. In order 
to understand the low-energy properties, it is natural to 
switch to fermionic language, performing a phase-shifted 
Jordan- Wigner transformation 



for i < M and 
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for i > M. Equation (23) then becomes: 
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which can, in turn, be rewritten as [see Fig. [3ja) 
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in terms of Majorana operators (A5). By (i) we denote 
here i ^ M and ... in Eq. (27) stand for terms of the 
form 7m7m+i, Jm7m+i, and jm7m+i, which involve 
the gapped states and thus have essentially no effect on 
the interaction of the end Majorana zero modes 7m and 
7m+i when t< 1. Making exchange interaction at the 
junction {M, M + l} anisotropic in the xy plane does not 
alter the essential outcome. For example, a spin exchange 
— t<r M <3" m+1 , where m' = (cos (//, sin<//, 0), results in the 
coupling ztcos <fi' cos(<p — (f>')jMjM+i m fermionic lan- 
guage, which is reduced by a geometric factor cos ci> f and 
phase-shifted by cj)' with respect to the isotropic XX case, 
Eq. (27). In particular, we see that the term <5"m^m+i 
in EqH23) corresponds to ci> f = tt/2 and thus does not 



couple Major anas 7m and 7m+i- 

According to Hamiltonian (27), the low-energy spec- 
trum is governed by the Majorana modes 71, 7 at, 7m, 
and 7m+i- The first pair of these defines a fermionic zero 
mode, b = (71 + £7at)/2, while the second pair fuses into 
an ordinary finite-energy fermion, a = (jm + Z7m+i)/2. 
The total Hamiltonian in the subspace of these four Ma- 
jor anas is given by 
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If the orientation <fi of the Ising coupling in the left spin- 
chain segment of our system is variable and controlled 
by some external handle (which, e.g., defines the shape 
of quantum dots that provide spin states), the latter ex- 
periences backaction torque 
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FIG. 3. (a) Two fused spin chains in the Majorana represen- 
tation. The end Majoranas 71 and 7/v form a fermionic zero 
mode b = (71 + z7at)/2, while the junction Majoranas 7m 
and 7M+1 fuse into an ordinary fermion a = (jm +^7m+i)/2, 
with energy e = 2tcos0. Here, t is the exchange coupling 
across the {M, M + 1} junction and cj) is the relative tilt of 
the Ising exchange interactions in the segments 1 < i < M 
and M + 1 < i < N. (b) When = 0, occupying fermion 
state a creates the spin domain wall sketched in the figure, 
while the domain wall corresponds to the empty fermion state 
when 6 — tt. 



in response to Hamiltonian (28). If the fermion state 
a is occupied, r z — Jtsin</>, while an empty state gives 
r z = —Jts'mcj). The ground- state configuration corre- 
sponding to Eq. (28) results in r z = — Jsgn(tcos (j)) sin (j). 



The period of the torque t z (4>) in equilibrium is thus 
half of that for a given (either empty or occupied) eigen- 
state. This is fully analogous to the fractional Joseph- 
son effect in topological semiconducting wiresP Consider 
the system initially (in the absence of coupling t) pre- 
pared in the state such that b\\i) = 62 K) = 0> where 
bi = (7i+*7m)/2 and b 2 = (7M+i+^7iv)/ 2 (i- e -> ordinary 
fermions composed of end Majoranas in the disjoint topo- 
logical segments). In the basis {1, a), b\ a^b^}\0), where 
a\0) = b\0) = 0, it is then given by \i) = (l+iattf)\0)/y/2. 
Turning t on then in Eq. ( [28] ) splits the degeneracy be- 
tween states |0) and a^^O), resulting in opposite torques, 
=F Jts'mcj) (which is zero on average in 

Torque measurement of two fused qubits (each corre- 
sponding to a topological segment of the spin chain) is 
thus sensitive to the initial two-qubit state, providing 
a possible read-out mechanism. While this is techno- 
logically challenging, recent strides in quantum measure- 
ments of nanomechanical cantilever d 22 * 23 * appear encour- 
aging to this end. In a flux qubit realization of the Ising 
chain, 24 on the other hand, such effective torque corre- 
sponds to a purely electronic signal. 

The most straightforward way to visualize fused qubits 
[see Fig. [3^b)] is to consider cases (j) = or 7r, which, ac- 
cording to Eq. (23), corresponds to collinear x-Ising seg- 

at the 

junction. (As noted above, the term —Jtcj y M c7 v M+1 



ments that are weakly coupled by — ^t^M^M+i 

(As noted above, the term —Jtcj v M cj y M+1 has 
no effect on the low-energy properties, to the leading or- 
der in t.) The ordinary fermion formed in the junction, 



Eq. (28), thus corresponds to an Ising domain wall with 



energy 2 Jt. The absence or presence of this domain wall 
can be revealed by simultaneous measurement of spin at 
the outside ends, i = 1 and i = N. In the presence 
of some Heisenberg interaction, furthermore, the domain 
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wall can move in response to an applied magnetic field 
along the chairP^J (inducing emf around the wire by 
the Faraday effect) pi which offers alternative schemes 
for information read-out as well as transmission. Finally, 
we note that well-established spin-to-charge read-outs^ 
can be applied to individual spin sites (such as quantum 
dots) constituting a spin-chain qubit, providing a simple 
mechanism for single-qubit measurements. 
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Appendix A: Generalities 



V. DISCUSSION 

External magnetic field fluctuations can in principle 
cause classical (i.e., flip) as well as phase errors in our 
logic space. Perhaps most problematically, the ferromag- 
netic order parameter (a x ) can couple to external mag- 
netic fields along the x direction, dephasing a single qubit 
in the {|=>), |^=)} basis or flipping it in the {|+),| — )} 
basisP This issue can be mitigated by adjusting the ex- 
change Hamiltonian of our spin chains in accordance with 
the structure of the dominant noise source, thus estab- 
lishing the qubits in a "decoherence-free subspace."^ As 
the simplest example, in the presence of long- wavelength 
magnetic-field fluctuations, encoding qubits with antifer- 
romagnetic rather than ferromagnetic spin chains should 
significantly reduce decoherence. Temperature, further- 
more, has to be much less than the quasiparticle exci- 
tation gap scale J. Despite mathematical similarity of 
topological spin-based qubits to their superconducting 
Majorana cousins, the dominant quantum errors there 
are of a very different physical origin.^ 



Here, we review introductory material that is neces- 
sary for constructing spin-chain based topological qubits. 
Consider the following general anisotropic ID spin- 1/2 
chain of N sites with nearest-neighbor exchange cou- 
plings and open ends: 
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The Hamiltonian H is normalized by the x-Ising ex- 
change constant J, and a, /?, h are dimensionless vari- 
ables parametrizing the relative strength of the remain- 
ing terms. In particular, if a = j3 = 0, we have the 
Ising chain in transverse field h. a = j3 — 1 describes 
the Heisenberg chain and a = 1, /3 = defines the XX 
model. The desirable anisotropic spin Hamiltonian can 
be accomplished utilizing spin-orbit interactions in elec- 
trostatically coupled lateral quantum dot d 3Q * 31 lor Wigner 
crystals.^ Effective spin chains can, furthermore, be en- 
gineered by an array of coupled superconducting flux 
qubits!^ 

Performing the Jordan- Wigner transformation^ 
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where af = (a x db ia y )/2 and hi = (af + l)/2 
defines fermionic representation for the Hamiltonian: 



c]ci 



H 



N-1 N-1 N 

= - E [( c l - c »)( c l+i + c *+i) - a (4 + a)(4+i - ci+i)] - p E - x )( 2 ^+i - x ) - ft E( 2 ^ - x ) • ( A3 ) 



2=1 



One easily verifies that {q,c^} = 5ij and {ci : Cj} = 
0, which is accomplished by the Jordan- Wigner string 
e ^E i <^i i ^ n otic e that the XX model gives for the 
first term in Eq. (A3) a simple hopping Hamiltonian 

N-1 

- 2 c U+i + H.c. , (A4) 

i=l 

while the Heisenberg model in addition has nearest- 
neighbor interaction ex r^n^ + i. The field term ex h in the 
fermionic language enters as a chemical potential shift. 
According to the form of the fermionic Hamiltonian 



(A3), it is convenient to define self-conjugate Majorana 



operatorPEUMl 



Ji = c\+Ci, ji=i(4-Ci), s.t. {7/,7j} = 2Jjj , (A5) 



which can be viewed as the real and imaginary parts of 
the fermionic fie ld op erators, respectively. The first term 
of Hamiltonian ( A3 ) then becomes 



N-1 



i ^{lili+i -«7i7i+i)- 



(A6) 
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1. Pure Ising chain 

In the special case of a pure Ising model, a = f3 = 
h = and thus H/J = i Y^iLi 1 7i7i+i? which means that 
the operators 71 and 7^ dropped out of the Hamilto- 
nian completely. Combining them into a new fermionic 
operator 



7i + ijN 



s.t. {M f } = 1 



(A7) 



and also pairing the remaining 2 (TV — 1) Majorana oper- 
ators 7's into (N — 1) fermionic operators = (7^ + 
z7^ + i)/2, such that H = a) Aa, with some Hermitian 
(TV — 1) x (N — 1) matrix A, we conclude that all eigen- 
states corresponding to A acquire double degeneracy 
due to b. Namely, for all eigenstates |0) a such that 
H\0) a = e a |0) a and 6|0) a = 0, there is also an orthogonal 
degenerate eigenstate |l) a = M|0) a . Most importantly 
for us, this concerns also the degenerate ground states 
|0) and |1). 

The end Majoranas that dropped out of the Hamilto- 
nian, in the case of a pure Ising model, give rise to double 
degeneracy of the ground state. This is of course obvious 
in the original spin representation, according to which we 
trivially have two ground states: 



->>) and |<=) 



where 



-f 1 

V2 \ 1 



— — , and \<-) = — = 



1 

3 



-1 



") , (A8) 



(A9) 



are the single-spin eigenstates of a x , and we assumed 
J > 0, to be specific. 

In addition to the time-reversal symmetry, the ground 
states (A8) spontaneously break the following Z2 sym- 
metry: 



N 



(A10) 



2=1 



which flips |=>) into |^=) and vice versa (flipping the asso- 
ciated magnetic order parameter (cr x )). In the fermionic 
language, 



N 



V = n(2n 4 - 1) = (-l)S«=i(« 



-1) 



(All) 



is the parity operator, with eigenstates ±1 corresponding 
respectively to the even/odd number of fermionic holes 
present in the chain. Since the ground states |0) and 
|1) have well-defined parity, they must correspond to the 
symmetric and antisymmetric combinations of |=>) and 
|<<=). Specifically, 



(A12) 



if N is even and vice versa if N is odd. In the spin 
language, the Majorana fermion operator b (such that 
6|1> = |0) and b\0) = 0) is 



(A13) 



where Z = n^i( — a i) 18 ^ ne s tring operator. Note that 
since b is a nonlocal operator because of the string Z, 
the degenerate ground states |0) and |1) cannot be inter- 
preted in terms of local spin states. The Z2 symmetry 
(A10), furthermore, stabilizes the degeneracy against lo- 
cal weak parity-conserving perturbations P 

It is important to point out a key difference between 
the aforementioned degeneracy in a spin chain vs a topo- 
logical superconducting wireP While the systems be- 
come mathematically equivalent in ID with open bound- 
ary conditions, it is not so once we switch to periodic 
boundary conditions.^ In particular, the Ising chain pre- 
serves degeneracy (which, according to the time-reversal 
symmetry, is unaffected by additional Heisenberg and/or 
Dzyaloshinsky-Moriya interactions, as long as some mag- 
netic ordering remains) , while the periodic boundary con- 
ditions in general lift the degeneracy of Majorana end 
states in the superconducting wires. This means that 
there is nothing fundamentally special about the end re- 
gions in our spin chains, and the Majorana-fermion anal- 
ogy with topological superconducting wires is only ful- 
filled for open boundary conditions. 

While both open and cyclic spin chains allow for the 
ground-state degeneracy, open chains are expected to be 
more robust to higher-order global |=>) O |^=) flipping 
processes than their closed counterparts. A closed chain 
can be flipped by a virtual local single-spin flip event 
(at an arbitrary location) that costs energy ~ 2 J, which 
then produces two domain walls that propagate around 
the chain and annihilate each other (thus returning the 
extra energy) after flipping the entire chain. In an open 
spin chain, on the other hand, the domain walls would 
need to be nucleated either at the ends or exactly at the 
center of the chain, or else the extra energy 2J would 
be returned at the two ends (one half of it at a time) 
at different times, thus requiring an additional virtual 
process. This observation reveals a topological character 
of our spin-chain logic-state subspace, which is distinct 
from the topological nature of the Majorana end states 
in superconducting wires. 



2. Ising chain in a transverse magnetic field 

Let us look more closely at the problem of Ising chain 
in a transverse magnetic field: 



H 

1 



N-l 



N 



(A14) 



While the time-reversal symmetry is broken, V still gen- 
erates a Z2 symmetry of Hamiltonian ( A14). In the ther- 



modynamic limit, the ground state remains degenerate 



8 



for | ft | < 1, while a paramagnetic phase with a unique 
ground state sets in for |ft| > 1. |ft| = 1 is thus the criti- 
cal point for a 2nd-order quantum phase transition.^ 



Diagonalizing Hamiltonian ( A14) is most convenient in 
the fermionic representation (A3), 



H 



N-l 



N 



noninter acting or, in other words, they have a vanish- 
ing correlation length £. A small finite magnetic field 
(and/or more general exchange interactions) gives rise to 
a finite length £, while still preserving the topological na- 
ture of the ground-state doublet of the wire whose length 
L ^> £. 16 Interaction between the end Majoranas is then 
generally given by 



1 = 1 



(A19) 



which is easily accomplished by the Bogolyubov trans- 
formation for the momentum eigenstatesP^ 



-ikj 



c k = -y= 2_j e C 3 ' ^ k cos VkCk + i sin VkC- k >. 

3 



where k = 2ir(n/N) G (— 7r,7r] and 



9 = - tan" 1 sinfe 

2 cos k + ft ' 



(A16) 



(A17) 



where t ~ e~^l L is the tunneling amplitude for a domain 
wall to propagate along the length of the chain, zipping 
\=>) into l^=) and vice versaP In terms of the complex 
fermions JA7L this Hamiltonian is given by 2t6 t 6, such 
that the lowest-energy splitting is 2t (while the other ex- 
citations are gapped approximately by the domain-wall 
energy of 2 J). 

Note that the sign of J (i.e., the ferromagnetic, J > 0, 
vs antiferromagnetic, J < 0, case) is of no consequence, 
since it can be flipped in Hamiltonian ( A14) by a unitary 



Here, for simplicity, we are impo sing periodic boundary 
conditions for Hamiltonian ( |A15| ), which does not affect 
bulk properties.^ In terms of these Bogolyubov quasi- 
particles, the Hamiltonian (in the thermodynamic limit) 
becomes (dropping an unimportant overall constant) 36 



transformation that amounts to a 7r rotation around the 
z axis (thus a z — >> a z and a x,y — )> — a x,y ) for all odd 
sites 

Finally, performing a global <fi rotation around the z 
axis, we obtain the following Hamiltonian: 



H 

7 



N-l 



N 



°i (J i+l- h Z^ (J i ' 



(A20) 



H = ^e/c7fc7/c , e fe = 2J\/l + ft 2 + 2ft cos k . (A18) where a z m = a L ■ m and m = (cos <j>, sin 0, 0). In the 
k fermionic language, Eq. (A2), this rotation generates a 

global gauge transformation: q — » Qe*^. The Ising in- 
teraction thus becomes 



The gap in the excitation spectrum closes at |ft| = 1, at 
which point we have a phase transition from the ordered 
phase with degenerate ground state (for |ft| < l), 16 to an 
ordinary quantum paramagnet phase with a nondegen- 
erate ground state (for |ft| > 1). Introducing exchange 
interactions along the z axis leads to nearest-neighbor 
interactions in the fermionic language, Eq. (A3). This 
needs to be avoided as strong interactions can sup- 
press the superconducting gap that stabilizes ground- 
state degeneracy.^ 

In the case of a vanishing field ft, the degeneracy of 
the ground state is exact, irrespective of the length of 
the chain. This means that the end spin Majoranas are 



'i+l 



(4e- 



•i+i 



C i C i+l e 



-2i<f> | 



H.c. 



) 

(A21) 



so that a <j) rotation in spin space is equivalent to the 2</> 
phase change of the superconducting order parameter, 
A — ^ Ae -2 ^, in the fermionic language. In particular, 
rotating the Ising axis by ir/2 [see Fig.[2|b)] is equivalent 
to changing the sign of the superconducting gap, which 
explains why the geometric manipulations in Fig. |2|b) 
mimic Majorana braiding (and the associated phase gate) 
in the superconducting wiresP 
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